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In  permafrost  areas,  in  addition  to  water  impermeable  dams,  use  is 
made  of  water-permeable  ones.  The  thermal  conditions  of  these  dams  have 
been  fairly  widely  studied  under  the  conditions  of  the  planar  problem. 


On  the  basis  of  the  functions  available  for  determining  the 
temperature  steady  state  of  water-permeable  dams,  this  paper  proposes 
dependencies,  utilizing  N.  K.  Girinskiy's  hydrodynamic  functions  [4]: 


Specific  flow  rate  potential 


<p  m3/day 


and  current  function 


m'Vday 


and  it  compares  the  calculation  results  obtained  by  a known  method  [3] 
with  the  results  of  calculations  performed  in  accordance  with  the 
proposed  functions. 

IVe  examine  a dam  made  of  homogeneous  isotropic  ground  sitting  on 
a layer  of  the  same  type  of  ground  as  in  the  body  of  the  dam,  which  is 
underlaid  by  a water  confining  stratum.  The  water  temperature  in  the 
reservoir  and  its  level  are  constant,  but  in  the  lower  race  there 
is  no  water. 


Girinskiy's  filtration  potential  for  an  unconfined  flow  in  homo- 
geneous ground  has  the  form  [1,  4]: 

h2  z 

y - ky-  in  / day  (1) 

where:  k m/day  - filtration  factor. 


h m - depth  of  filtration  flow. 


We  use  Girinskiy's  potential  to  plot  a depression  curve.  For  this 
purpose  potentials  sj  and  ^ corresponding  to  the  depths  of  the 

filtration  flow  in  the  upper  race  h^  and  at  the  drain  level  h^  are 


designated  in  relative  units  as  sP*  =1  and  = 0;  then  the  depth 

lg  2g 

of  the  flow  in  vertical  lines  can  be  obtained  by  the  following  formula, 
provided  that  the  distance  between  the  cross-sections  where  ^ and 

were  determined  is  equally  divided: 


/i  - h.  TUT-  ~fif) 


+ n. 


(2) 


where  * * - the  relative  potential  of  the  intermediate  vertical  line  in 
fractions  of  the  difference  between  the  upper  and  lower  races. 


A depression  curve  plotted  according  to  the  values  of  (2)  is  a 
Dupuis  parabola. 

The  filtration  speeds  at  any  point  on  the  vertical  line  corresponding 
to  the  flow  depth  are  equal  and  the  vertical  component  is  equal  to  zero 
[4]. 


The  Girinskiy  potential  is  used  to  calculate  the  specific  flow  rate 
of  the  filtration  flow  [2]  . 


g 


m^/day 


(3) 


where  is  the  distance  between  the  cross-sections  where  $ 
are  determined. 


lg 


and  v 


2g 


The  specific  flow  rate  makes  it  possible  to  calculate  the  vertical 
speeds  from  the  function: 


v = j*-  m/day.  (4) 

The  vertical  flow  is  divided  by  surfaces  which  are  close  together 
and  horizontal  in  p layers,  and  then  the  height  of  each  layer  is 


A A =■  ~s~  ■ u (5) 

In  expressing  filtration  flow  in  Girinskiy's  hydrodynamic  functions, 
the  potential  increment  for  two  adjacent  vertical  lines  separated  by  a 
distance  A Vt  can  be  represented  in  the  following  form: 

c 6-  & m3/day-  (6) 

Heat  exchange  in  water-permeable  ground  is  determined  for  steady 
state  conditions  by  the  Fourier-Kirchhof  differential  equation  [3], 
which  by  means  of  the  Girinskiy  potential  for  vertical  cross-sections 
of  flow  can  be  presented  in  the  form: 


(7) 


At  considerable  filtration  speeds  when  the  convective  transfer  of 
heat  along  the  current  line  predominates  over  conductive  transfer,  we 

7 

ad 

can  ignore  the  term  ha^  - — 0 in  expression  (7) , and  then  the  Fourier- 

X '‘’g*' 

-Kirchhof  equation  has  the  form: 


In  this  form  the  equation  is  solved  numerically  in  finite  differences. 
In  the  vertical  flow  a point  called  the  central  is  selected;  the  tempera- 
ture of  this  point  is  designated  as  0 ; and  <9  jQw  indicate  the 

temperatures  of  the  points  in  front  of  and  behind  the  central  in  the 
direction  of  the  filtration  flow;  ^ houn(j  and  ^ bound  s*10w  t*ie 

temperatures  of  points  above  and  below  the  central  on  one  vertical  line. 

In  equation  (8)  the  derivatives  in  the  finite  differences  are 
presented  in  the  following  form: 


By  substituting  these  values  into  (8)  with  allowance  for  (4),  (5) 
and  (6),  we  obtain  the  formula  for  computing  the  central  point  tempera- 
tures : 


to  (10),  the  temperatures  are  determined  at  the  points  which  divide  the 
flow  depth  into  the  accepted  number  of  layers. 


A determination  of  a dam's  temperature  state  according  to  (10)  is 
examined  using  the  example  of  designing  a dam  made  of  homogeneous 
isotropic  ground  with  a filtration  factor  of  K = 1.7  m/day  underlaid  by 
a horizontal  water  confining  stratum  whose  temperature  is  0°C. 

Above  the  stratum  is  a layer  of  the  same  ground  as  is  in  the  dam  body; 
it  is  5 m thick.  The  dam  has  tubular  drainage,  the  upper  slope  of  the 
dam  is  1:2  and  the  lower  is  1:1.5.  The  water  temperature  in  the  reservoir 
is  a constant  3°C,  the  water  depth  in  the  upper  race  is  h = 10  m,  in  the 
lower  race  there  is  no  water. 

The  Girinskiy  potential  is  used  to  plot  a depression  curve:  in  the 
flow  10  vertical  lines  4.5m  apart  are  assumed  and  the  speeds  on  the 
verticals  are  determined  from  (4).  The  results  are  summarized  in  the 
table  (Table  1).  bach  of  the  vertical  lines  was  divided  into  three 
equal  sections,  and  temperatures  were  determined  at  the  division  points 
and  at  the  extreme  upper  points  of  the  vertical  lines;  this  corresponds 
to  the  temperature  state  on  the  depression  surface.  In  this  case,  upor 

extension  of  each  vertical  line  at  distance  Ah  = y from  the  depression 

surface  the  temperature  was  assumed  to  be  zero,  and  therefore  the 
temperatures  at  points  at  distance  Ah  from  the  stratum  surface  and  on 
the  depression  surface  are  symmetrical . 


Isotherms  have  been  plotted  from  the  calculated  temperatures 
(Figure  1).  This  same  figure  indicates  the  temperature  values  at  the 
comers  of  the  curvilinear  quadrants  of  the  hydrodynamic  grid,  as 
calculated  from  (11). 


A comparison  indicates  that  the  temperature  values  at  the  corresponding 
points  of  the  dam  cross-section  coincide  well;  this  shows  that  the  method 
of  approximately  calculating  filtration  through  a dam  can  be  applied  to 
determine  the  temperature  steady  state  of  water-permeable  dams. 

The  calculated  temperatures  in  the  dam  cross-section  are  used  to 
determine  the  temperature  gradients  from  the  thawed  ground,  and  the 
thawing  of  the  frozen  base  during  the  first  five  years  is  determined 
with  the  aid  of  the  Stefan  equation. 

Since  the  frozen  ground  has  a temperature  of  zero,  the  gradients 
from  the  frozen  ground  are  equal  to  zero;  in  this  case  thawing  is  deter- 
mined as : 


AV.g“  ldeg- 


where : 


M cal 


X — j j — is  the  heat  conductivity  factor  of  the  thawed  ground; 

tg  m*deg*day  J & ’ 

id  $ deg  is  the  modulus  of  the  temperature  gradient  from  the 
| 8 tg  m thawed  ground; 

y^  n L — is  the  latent  heat  of  the  change  in  aggregate 

m moisture  state  in  the  ground  pores; 


At  day 


is  time. 


When  Xtg  = 0.0544  ^day-deg  and  Yf  n L = 27 -8  the  intensity  of 

thawing  for  a period  of  five  years  is  determined  (the  values  of  A 
are  given  in  the  table) . 

TABLE  1 . 


No.  of  vertical 



lines 

Flow  depth 

u h 

1 15,0 

1 

t 

14,3 

- 


▼ 


Conclusion. 

In  making  temperature  calculations  of  water- permeable  dams,  in 
addition  to  precise  filtration  computation  methods  it  is  also  possible 
to  use  an  approximate  method  of  calculation  which  utilizes  Girinskiy's 
hydrodynamic  functions. 
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